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Abstract. I introduce two family of exactly solvable models for multiatomic hetero- 
nuclear and homo-nuclear molecular Bose-Einstein condensates through the algebraic 
Bcthe ansatz method. The conserved quantities of the respective models are also 
showed. 



1. Introduction 

One of the most interesting recent experimental achievements in physics is the one that 
led to realizations of Bose-Einstein condensates (BEC), by taking dilute alkali gases to 
extremely low temperatures [IJ[2]. Since then, a great effort has been devoted to the 
comprehension of new phenomena involving this state of matter as well as its properties, 
either experimentally or theoretically. On the experimental side I could mention a 
molecular BEC compound that has been obtained combining different techniques [3], 
leading this kind of research also in the direction of a chemistry of BEC, where, for 
instance, by Feshbach resonances [IH6] or photo-association [7J[8] the atomic constituents 
may form molecules. 

Many compounds of diatomic homo-nuclear molecular BECs [9HTT] have been 
produced since the first realization . Also, diatomic hetero- nuclear molecular 
BECs have been detected using these techniques [T3HT9] . Actually, due to the rapid 
technological developments in the field of ultra-cold systems, it is believed that some of 
these experiments may be just the dawn of the study of multiatomic molecules [201 121] . 
More recently the experimental evidence for Efimov states in an ultra-cold cesium 
gas [22H21] and a mixture of ultra-cold potassium and rubidium gases [25] provides 
a physical ground for the investigation of triatomic and tetratomic homo-nuclear and 
triatomic hetero-nuclear molecular BECs. 

These results boosted the search for solvable models that could describe some of the 
BEC properties [26T - H2] . The rationale beneath these studies is that through exactly 
solvable models it is possible to fully take into account quantum fluctuations, going 
beyond the usual mean field approximations. Then, I expect that this approach may 
provide some impact in this well as a contribution to the field of integrable 
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systems itself [121111] • In this paper I will use the algebraic Bethe ansatz method. The 
algebraic formulation of the Bethe ansatz, and the associated quantum inverse scattering 
method (QISM), was primarily developed by the group of mathematical physicists in 
St. Petersburg [15T - H9] . The QISM could be used to study the one-dimensional spin 
chains, quantum field theory in one-dimensional bosons interacting systems [50] and 
two-dimensional lattice models [51], systems of strongly correlated electrons [52|l53]. 
conformal field theory [31], as well as precipitated the notion of quantum algebras 
(deformations of universal enveloping algebras of Lie algebras) [55ti58] . For a pedagogical 
and historical review see [39] . 

Owing to recent insights in the understanding of the construction of Lax operators 
it is possible to obtain solvable models suitable for the effective description of the 
interconversion interactions occurring in the BEC. Inspired by some of these ideas I 
present, in the present paper, the construction of two complete family of Bethe- ansatz 
solvable models for both homo-nuclear and hetero-nuclear molecular BECs obtained 
through a combination of three Lax operators constructed using special realizations of 
the su(2) Lie algebra and of the Heisenberg-Weyl Lie algebra, as well as a multibosonic 
representation of the si (2) Lie algebra, discussed recently in [60]. Notice that the models 
obtained through this construction do not have spatial degrees of freedom. 

The paper is organized as follows: In Section 2, I will review shortly the algebraic 
Bethe ansatz method and present the Lax operators and the transfer matrix for both 
models. In Section 3, I will present a family of multiatomic homo-nuclear models and 
their solutions. In Section 4, 1 will present a family of multiatomic hetero-nuclear models 
and their solutions. In Section 5, I will make my remarks. 

2. Algebraic Bethe ansatz method 

In this section we will shortly review the algebraic Bethe ansatz method and present 
the transfer matrix used to get the solution of the models [28,61]. We begin with the 
g/(2)-invariant i?-matrix, depending on the spectral parameter u, 



\ 1 / 

with b(u) = u/(u + if) and c(u) = i]/(u + 77). Above, r\ is an arbitrary parameter, to be 
chosen later. It is easy to check that R(u) satisfies the Yang-Baxter equation 



Here Rjk{u) denotes the matrix acting non-trivially on the j-th and the k-th spaces and 
as the identity on the remaining space. 

Next we define the monodromy matrix T(u), 



R(u) 



( 1 \ 

b(u) c(u) 
c(u) b{u) 



(1) 



R 12 (u - v)R 13 (u)R 23 (v) = R 23 {v)Ri 3 {u)R l2 {u - v). 



(2) 




(3) 
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that satisfy the Yang-Baxter algebra, 

R 12 (u - ujTiC^TzCv) = T 2 (v)T 1 (u)R 12 (u - v). (4) 

In what follows we will choose different realizations for the monodromy matrix 
7t(T(m)) = L(u) to obtain solutions of two family of models for multiatomic hetero- 
nuclear and homo-nuclear molecular BECs. In this construction, the Lax operators 
L(u) have to satisfy the relation 

Ri 2 (u - v)L 1 (u)L 2 (v) = L 2 {v)Li(u)Ri 2 {u - v). (5) 

where we use the notation, 

L 1 = L{u)®I and L 2 = I®L{u). (6) 
Then, defining the transfer matrix, as usual, through 

t(u) = tm(T(u)) =7r(A(u) + D(u)), (7) 

it follows from (j4j) that the transfer matrix commutes for different values of the spectral 
parameter; i. e., 

[t(u),t(v)} = 0, V«, v. (8) 

Consequently, the models derived from this transfer matrix will be integrable. Another 
consequence is that the coefficients Ck in the transfer matrix t(u), 

t(u) = Y,CkU k , (9) 

k 

are conserved quantities or simply c-numbers, with 

[Cj,C k ] = 0, Vj,k. (10) 

If the transfer matrix t(u) is a polynomial function in u, with k > 0, it is easy to 
see that, 

1 d k t(u) 



C = t(0) and C k 



(11) 

u=0 



k\ du k 

We will use three solutions of the equation (JHJ): 
(i) The L s [u) Lax operator: 

in terms of the su{2) Lie algebra with generators S z and S^ 1 subject to the 
commutation relations 

[S z , S+] = ±S ± , [S + , S-] = 2S Z . (13) 
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(ii) The L?(u) Lax operator: 

U(«)=[ U+ P /:)• (14) 

in terms of the Heisenberg-Weyl Lie algebra with generators Nj, j, j' and /, subject 
to the commutation relations 

[NjJ] = -j, = and [/,*] = 0, (15) 

where * means Nj, j or j* . 

(iii) The L A {u) Lax operator: 

in terms of the si (2) Lie algebra with generators A and A ± , subject to the 
commutation relations 

[A^,A + ] = A , [A ,A ± ]=±2A ± . (17) 

Using the co-multiplication properties of the Lax operator and the gl{2) invariance 
of the i?-matrix, we can obtain different realizations for the monodromy matrix: 

(i) Multiatomic homo-nuclear: For the multiatomic homo-nuclear molecular BEC 
model we choose 

7r(2») = L{u) = r]- l GL j (u -5- r ] ~ 1 )L A {u + uj), (18) 
with, G = diag(—, +), from which we find the following transfer matrix 



V 

+ r]~ 2 (u + u-^A ) + bA + + b^A^, (19) 
with 

t(0) = rf x u{b + 2rf x ) - uN h + i(5 - r]N b )A + bA + + (20) 
and, discarding c-number terms, the conserved quantities are, 



t(u) = -r]- 1 (u + u+ +A )(u - 5 - r]- 1 + r]N b ) 



Co = t(0), 

Ci = + iV 6 . 



(21) 
(22) 
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(ii) Multiatomic hetero-nuclear: For the multiatomic hetero-nuclear molecular 
BEC models we choose 

tt(T(«)) = L(u) = r)- 1 u-GL s (u~)L A (u + ), (23) 

with = u±u, G = diag(+, — ), from which we find the following transfer matrix 

t(u) = u~A - 2u + S z + ri(S + A + + S'AJ), (24) 
and the conserved quantities, 

Co = t(0), (25) 
d = Ao- 2S Z . (26) 

In the next sections we will describe the models and its integrability by the algebraic 
Bethe ansatz method, using different realizations of the algebras (fl3"|) . ffl~5l) and (TTT|) . 
The Hamiltonians are written in the Fock space using the standard notation. We are 
considering the coupling parameters real, such that the Hamiltonians are Hermitian. 
In the diagonal part of the Hamiltonians, the Uj parameters describe the atom-atom, 
atom-molecule and molecule-molecule S'-wave scatterings and the fij parameters are the 
externals potentials. The operators Nj are the number operators of atoms or molecules. 
In the off diagonal part of the Hamiltonians the parameter Q is the amplitude for 
interconversion of atoms and molecules. 

3. Multiatomic homo-nuclear molecular models 

In this section we present the integrability of a new family of Hamiltonians describing 
multiatomic homo-nuclear molecular BEC. The Hamiltonians that describes the 
interconversion of homogeneous molecules labelled by b with / atoms of type a are 
given by 

H = U a N 2 a + U b N% + U ab N a N b + n a N a + fi b N b 

+ n((a t ) J b a-(N a ) + a-(N a ) b ] (a) 1 ), (27) 

where a_(N a ) is a function of N a that controls the amplitude of interconversion Q. 
This indicates that the density of atoms N a has some influence in the generation of a 
bound-state composed by I identical atoms. The I = 3 case was studied in (35]. The 
total number of particles N = N a + lN b is a conserved quantity. 
There is a multibosonic realization of the si (2) Lie algebra [60] 

A = a (N) } A- = a-(N)a l , A + = (a^) l a4N), (28) 



with 
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a (N) =j(N-R) + a (R), (29) 



AM 1 1 
"- ( - V) = \l^ TJ y( 1 (N-R) + a (R))( 1 (N-R) + l), (30) 

where N = a)a and I G N. The operator R is 

for 1 = 1, 

R= { l -l i± e -(^m/l)N (31) 

h L^woTT for Z>1 > 

m=l 

and acts on the states {\n)} as R\n) = n mod l\n). The function ceo(R) is a positive 
function of the spectrum of R defined by initial conditions. For n = r < I, we have 

y(JV-i2)|r) = 0|r), (32) 

with A = a>o(R) such that ao(R)\r) = Qo(r)\r) and R\r) = r\r). 

Now we will use this realization to show how to construct the Hamiltonian (1271) 
from the transfer matrix ffl9l and present their exact Bethe ansatz solution. It is 
straightforward to check that the Hamiltonian ( 1271) is related with the transfer matrix 
t(Q) (1201) . or with the conserved quantity C (12TT) . through 

H = n t(0), (33) 

where we have the following identification 



l 2 Ua ~ lUab + Ub 

n 



(34) 



9 = 2lU ab - 4U b , £ = 2/ 2 /ia - 2lfi b , (35) 

2lQ(u + 6) = (20?7 + ^)A^ - Vlpr] + £, (36) 

Al 2 nu57] + 8l 2 ttu = 7] 2 [n P 2 7] + AU b N 2 - 6Np 

+ 4J(£k) + ii b )N + (2lSk;-£)p], (37) 

with p = p(R) = la {R) - 2R. 

It is easy to see that p is a conserved quantity using the total number of atoms, N, 
to write the conserved quantity C\ f l2~2"|) as, 

C^jiV+lp. (38) 
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We can apply the algebraic Bethe ansatz method, using as the pseudo- vacuum the 
product state (|0) = |0)& <S> |r)^, with \0) b denoting the Fock vacuum state and \t)a 
denoting the lowest weight state of the s/(2) Lie algebra, where r = 0, 1, ... , / — 1, are 
the eigenvalues of R for iV = nl + r, with n G N, to find the Bethe ansatz equations 
(BAE) 

(l-^ + ^ + .lla,(r)) = fj m i,j = l,...,M, (39) 
Vi + u-ja {r) f^Tj v i ~ v j + V 

and the eigenvalues of the Hamiltonian ( 1271) . 

M 



E = n v - 1 (5 + V - 1 ) (w + |oo(r))lJ 



^- 2 (--^ (r))n^- (40) 



2 v V J-- 1 - w 

i=l 

The parameters 5 and cu are arbitrary and can be chosen conveniently. In the limit 
without scatterings, Uj — > 0, the BAE (l39i) can be write as, 

> = — rr- > --5, (41) 

+ u ao(r) f-f "o r V ; 

and for u = the eigenvalues (I4TJ1) become, 



£ = ( -a (r) + MJ JW- fijjui. (42) 

Now, the relation between the interconversion parameter and the externals potentials is 
simply, 

n = ^V^- (43) 

o 

4. Multiatomic hetero-nuclear molecular models 

In this section we present the integrability of a new family of Hamiltonians describing 
multiatomic hetero-nuclear molecular BEC. The Hamiltonians that describes the 
interconversion of heterogeneous molecules labelled by c with I atoms of type a and 
one atom of type b are given by 



H = U a N 2 a + U b Nl + U C N 2 C + U ab N a N b + U ac N a N c + U bc N b N c 
+ fi a N a + p, b N b + fi c N c 

+ n((o)) 1 6 f c a_(JV ) + a_(A" a ) c f b (a) 1 ), (44) 
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where a_(N a ) is a function of N a that controls the amplitude of interconversion Q. In 
the same way of the Hamiltonians ff271) . this indicates that the density of atoms N a has 
some influence in the generation of a bound-state composed by / identical atoms. 
The imbalance between the number of atoms a and the number of atoms b, 

J ab = N a - lN b , (45) 

is a conserved quantity and the total number of atoms, N = N a + N b + (I + 1)N C , can 
be writes with the other two conserved quantities 

l ac = N a + IN C} (46) 

X bc =N b + N c . (47) 

Using J ab and N, the 5-wave diagonal part of the Hamiltonian (jSJ) can be writte 

as, 

aj 2 b + (3N 2 + jNJ ab , (48) 
where we have used the following identification for the coupling constants 

U a = a + (3 + ~f, U b = al 2 + (3-~fl, U c = (3(l + 1) 2 , (49) 

U ab = - 2la + 2(3 - 7 (Z - 1), U ac = 2(3(1 + 1) + 7 (Z + 1), (50) 

U be = 20(l + l)-jl(l + l). (51) 
Now, using the following realization for the su(2) Lie algebra, 

S+ = tfc, S~ = c\ S z = b c , (52) 
and the multibosonic realization of the si (2) Lie algebra ( 1281) . 

A = a (N) } A- = a-(N)a l , A+ = ((f) l a-(N), (53) 

with 

a (N) =j(N-R)+a (R), (54) 



TV' 1 1 
n - m J^—^ijiN-^ + aoiR^jiN-Rj + l), (55) 



where N = a) a and I £ N, it is straightforward to check that the Hamiltonian 
is related with the transfer matrix t{u) (|24"|) . or with the conserved quantity Co (|25|) if 
u — 0, through 

H = cr + aJ 2 b + (3N 2 + 1 NJ ab + t(u), (56) 
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where the following identification has been made for the parameters 



/i a = 2—, fi c = -fi b = u + , Q = rj, a = --j~P, (57) 

with p = p(R) = lato(R) - 2R. 

We also can use the conserved quantities J a b and X ac to write the conserved quantity 
d (EE]) as, 

Ci = ^{Jab+X ac ) + -p, (58) 

showing that p is also a conserved quantity. 

We can apply the algebraic Bethe ansatz method, using as the pseudo-vacuum the 
product state (|0) = \t)a ® \<fi)-> with \t)a denoting the lowest weight state of the sl{2) 
Lie algebra where r = 0, 1, ... , I — 1, are the eigenvalues of R for iV = nl + r, with 
n E N and \(f>) denoting the highest weight state of the su(2) Lie algebra with weight 
m z ), to find the Bethe ansatz equations (BAE) 

(,. w _ rt + w+Wr)) Zvi-Vi-V iJ = 1 ^ Mim 



(vi-u + r]m z )(vi + bj - ^ao(r)) Vi-vj + r]' 



and the eigenvalues of the Hamiltonian 



E = a + a J 2 ah + /3iV 2 + ^NJ C 



ab 



M 



+ (u-u- r]m z )(u + u + -a (r)) J | 



77 , . . T-r U — V i + 1] 



U-Vi 

1=1 

M 

7] T-r U — Vi — T] 



(„ _ u + + w _ i ao ( r )) rr — 1. (60) 



The eigenvalues f )60|) are independent of the spectral parameter u and of the 
parameter u, that are arbitrary. 



5. Summary 

I have introduced two new family of multiatomic molecular BEC models for homo- 
nuclear and hetero-nuclear molecules and derived the Bethe ansatz equations and the 
eigenvalues. The conserved quantities are also derived. The multiatomic homo-nuclear 
and hetero-nuclear molecular BEC models were obtained through a combination of Lax 
operators constructed using special realizations of the su(2) Lie algebra and Heisenberg- 
Weyl Lie algebra, as well as a multibosonic representation of the si (2) Lie algebra. The 
dependence of the parameters with the size of the molecules is explicit. 
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